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ABSTRACT 

The  stability  of  the  critical  state  for  flexural  buckling 
of  a  plate  of  incompressible  isotropic  elastic  material  with 
arbitrary  strain-energy  function  is  studied.  Most  of  the  analy¬ 
sis  is  carried  out  without  restriction  on  the  magnitude  of  the 
aspect  ratio  of  the  plate.  However,  the  final  result  is  limited 
to  the  case  when  the  aspect  ratio  is  such  that  terms  of  fourth 
degree  in  it  may  be  neglected  in  comparison  with  terms  of  zero 
degree. 


2. 


1.  Introduction 

In  the  present  paper  we  consider  a  plate  of  incompressible 
isotropic  elastic  material  to  be  situated  with  its  edges  parallel 
to  the  axes  of  a  rectangular  cartesian  coordinate  system  x  .  The 
plate  is  acted  on  by  dead- loads  applied  normally  to  the  faces  of 
the  plate  which  are  perpendicular  to  the  1  and  3- axes  of  the  system 
X  .  The  faces  perpendicular  to  the  2 -axis  are  force- free.  The 
constraints  on  the  1  and  3- faces  are  such  as  to  permit  the  plate 
to  undergo  pure  homogeneous  deformations.  It  is  supposed  that  the 
load  in  the  1-direction  is  a  thrust.  As  this  thrust  is  increased 
static  bifurcation  solutions  in  the  12 -plane,  superposed  on  a  uni¬ 
form  extension  in  the  3-direction,  become  possible  at  certain 
critical  values.  These  solutions  may  correspond  to  buckling  of  the 
plate  of  the  flexural  or  barreling  type.  In  [1,2]  the  compression 
ratios  in  the  1-direction  at  which  these  critical  values  of  the 
thrust  are  reached  was  calculated  for  an  arbitrary  strain-energy 
function. 

In  a  previous  paper  [3]  we  discussed  the  stability  of  the 
states  of  pure  homogeneous  deformation  at  which  these  bifurcations 
occur,  with  the  assumption  that  the  strain-energy  function  is  neo- 
Hookean.  The  stability  criterion  employed  and  the  procedure  adop¬ 
ted  was  essentially  that  due  to  Koite.  An  equilibrium  state 

is  regarded  as  stable  or  unstable  accurc^.ii^ly  as  the  potential 
energy  of  the  system,  consisting  of  the  body  and  loads,  has  a 
proper  minimum  at  this  state  with  respect  to  all  infinitesimal 
deformations  satisfying  the  kinematic  constraints.  It  was  seen  in 
[3]  that  a  state  of  pure  homogeneous  strain  for  which  a  bifurcation 
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solution  exists  is  one  of  neutral  equilibrium.  The  state  is 
stable  if  the  potential  energy  of  the  system  is  smaller  for  this 
state  than  it  is  for  every  state  in  its  neighborhood  which  satis¬ 
fies  the  kinematic  constraints;  otherwise  it  is  unstable. 

With  this  criterion  it  was  shown  in  [3]  that  at  critical  com¬ 
pression  ratios  for  flexural  buckling  the  homogeneous  state  is 
stable  provided  that  the  aspect  ratio  C^-dimension/l-dimension)  is 
less  than  about  0.2,  and  unstable  otherwise.  This  implies  that 
the  immediate  post-buckling  behavior  will  be  stable  for  the  lower 
aspect  ratios  and  of  the  snap- through  type  at  the  higher  aspect 
ratios.  At  critical  compression  ratios  for  buckling  of  the  barrel¬ 
ing  type,  the  homogeneous  state  is  stable  for  all  aspect  ratios. 
These  calculations  were  carried  out  in  [3]  for  aspect  ratios 
ranging  from  zero  to  infinity.  The  formula  on  which  these  calcu¬ 
lations  were  based  is  extremely  complicated  and  acquires  meaning 
only  as  the  result  of  the  numerical  computations.  Accordingly, 
for  the  case  of  flexural  buckling  an  asymptotic  calculation  valid 
for  small  values  of  the  aspect  ratio  was  also  carried  out.  As  the 
aspect  ratio  tends  to  zero  this  result  agrees  with  the  classical 
result  of  Euler  based  on  the  theory  of  the  elastica. 

In  the  present  paper  analogous  calculations  are  carried  out 
for  an  incompressible  isotropic  elastic  material  with  arbitrary, 
rather  than  neo-Hookean,  strain-energy  function.  While  much  of 
the  analysis  (up  to  §7)  is  carried  out  with  no  restrictions  on 
the  magnitude  of  the  aspect  ratio,  the  final  result,  which  is  given 
in  equation  (11.25),  is  obtained  only  for  the  case  of  small  aspect 
ratio. 


Statement  of  the  problem 

We  consider  a  rectangular  plate  of  incompressible  isotropic 
elastic  material,  which  has  its  edges  parallel  to  the  axes  of  a 
rectangular  cartesian  coordinate  system  x.  Let  C  be  the  vector 
position,  relative  to  the  origin  of  the  system  x,  of  a  generic 
particle  of  the  plate  in  its  undeformed  state  (state  0)  and  let 
its  bounding  surfaces  in  this  state  be  the  planes 

Ca  *  ±  (A*1,2,3)  .  (2.1) 

We  suppose  that  the  plate  is  maintained  in  an  equilibrium 
state  of  pure  homogeneous  deformation  (state  I),  with  extension 
ratios  principal  directions  parallel  to  the  co¬ 

ordinate  axes,  by  uniformly  distributed  normal  tractions  applied 
to  the  surfaces  “  *  ■^3»  surfaces 

Cg  =  ±  ^2  being  force-free.  Let  and  II23  be  the  tractions, 

measured  per  unit  undeformed  area,  applied  to  the  faces  =» 
and  ?3  =  respectively. 

We  assume  that  the  surfaces  initially  at  »  ±  and 
^3  ■  ±  -tj  are  constrained  so  that  they  move  parallel  to  the  1 
and  3-axes  respectively,  but  points  on  them  are  free  to  move  in 
the  planes  normal  to  these  directions  (i.e.  the  tangential  trac¬ 
tions  on  these  surfaces  are  assumed  to  be  zero) . 

Let  X  be  the  vector  position  in  state  I  of  the  particle 
which  has  vector  position  C  in  state  0. 


s. 


* 

Then, 

Xa  »  XaCa  (A=1,2,3),  X^A2X3  -  1  •  (2.2) 

Now  suppose  that  the  plate  undergoes  a  further  deformation 
which  consists  of  a  uniform  stretch  in  the  3 -direction  and  a 
plane  deformation  in  the  12 -plane.  We  call  the  resulting  state 
of  deformation  of  the  plate  state  II.  and  we  write 

X  »  X  +  u  ,  (2.3) 

where 


Ui  =  “i(?1»52^  »  “2  *  *  “3  “ 

and  E  is  a  constant. 

Since  the  material  is  incompressible,  det||xj^  I  ”  ^  • 
It  follows  with  C2.2)-C2.4)  that 


(1+E)CXj^U2  2*^2'^- 


Let 

matrices 


5  »  I  ’  MXijll  Finger  strain 

in  states  I  and  II  respectively.  Then,  with  (2.2), 


'AB 


^A^AB 


^i,m*j ,m 


(2.6) 


*  Throiighout  this  paper,  latln  subscripts  take  the  values  1,2,3  and  greek 
subscripts  take  the  values  1,2.  The  Einstein  summation  convention  applies 
to  lover  case  latln  and  greek  subscripts  but  not  to  upper  case  latln  sub¬ 
scripts.  Also,  the  operators  3/3Cjj»  ““e  denoted  ,a  and  ,m 

respectively. 
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We  introduce  the  notation 

X  “  £  *  S  »  s  "  M  •  (2.7) 

From  (2.4),  (2.6)  and  (2.7),  we  obtain 


^11 

*  (2Xi+Ui^l)Ui 

.1*" 

2  c  = 

1.2  »  22 

(^^2*^2,2^S,2‘^S,1  ’ 

^33 

=  X^E(2+E)  ,  c 

12“ 

Si  = 

.l)S,l^fS-^S,2)S,2' 

(2.8) 

c.  . 
11 

=  0  (ij«23,32 

,13. 

31)  . 

Let 

^1*'2 

invariants  of  - 

i  and  I^,  Ig  those  of 

C  defined  by 

ii  * 

=  tr  Y  =  »  ^2  *  ii(^^  Y^}  *  ^2*  ^2  • 


With  the  notation 


X  -  X^/X^  , 


(2.10) 


"1“’  ■  “<..6“o,e  *  ’ 


it  follows  from  (2.5)-(2.9)  that 


“  i.  “  ^3Ci+j)  . 


(2.11) 


where 
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(2.12) 


2{X^(1-x2)u3^^^  +  (x2-x2)E}  +  <[u] 

+  (x2+2x2)e2  -  2x|e^(1+E)-^  , 
E{2(x2+x2-x2-xJx^)  +  (2+E)i 
+  E(x2+x2-5x2+3X^X^) 

-  e2[x2(4+E)  +  xJx^(4+3E)(1+E)~2]} 


Let  W  and  w  denote  the  strain  energies  per  unit  volume 
in  states  I  and  II  respectively.  Then, 


w  =  ,  W  =  wCI^.Ig) 


(2.13) 


We  introduce  the  notation 


3w 

TT  ’ 


*  3T;5i“ 


(2.14) 


=  Wqj(I3_,I2)  ,  W^o  -  Wjjg(I^,l2) 


The  increase  in  the  strain-energy  of  the  plate  in  passing  from 
state  I  to  state  II  is  given  by* 


(w-W)dC,d?2  *  2^3  I  (W„i„4  W^gVp 


+  W^^^  +  +  wf5^)d53^d52  » 


(2.15) 


*  The  domain  of  integration  for  the  double  integrals  is  the 
rectangle  ^  . 
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where 


i  W  .  .i  i_i  i-  , 
24  agyo  a  B  y  o 


(2.16) 


and  denotes  the  usual  Taylor  series  remainder. 

Let  n  .  denote  the  Piola-Kirchhoff  stress  in  state  I  . 
ai 

Then  (cf.[l],  equation  (3.6)) 


"ll  “  2A^a-A2)(W^+x2j^2)  , 

^33  “  2CX3-x2a-1)(W^+x2w2)  ,  (2.17) 

=0  (ai  ^  11,33)  . 


It  follows  from  the  constraint  conditions  on  the  surfaces 
C,  =  that  the  displacement  field  u  must  satisfy  the 

conditions 


Ul^2(±^l,52)  =  0  .  (2.18) 

With  the  further  assumption  that  the  displacements  of  these  two 
surfaces  are  equal  and  opposite,  we  obtain 

»  (2.19) 

where  e  is  a  constant. 

The  resultant  forces  acting  on  the  faces 
“*^3  state  I  are  obtained  from  (2.17)  as  ±  and 

±  ^®sP®ctively .  We  suppose  that  state  II  is  reached 

from  state  I  under  dead-loading  conditions.  Then  these  are  also 

the  resultant  forces  acting  in  state  II.  The  increase  in  the 

potential  energy  of  these  forces  in  the  passage  from  state  I  to 

- .  . ^ ■  A 
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state  II  is 

•  C2.20) 

With  (2.17)  and  (2.19)  this  may  be  written  in  the  form 

-U3  I  |{A^(l-l2)(;w^+x2w2)u^^^+(A2-X^)(W^+A2w2)E}dC^d52.  (2.21) 

The  increase  G[u]  in  the  total  potential  energy  of  the 
system,  consisting  of  the  body  and  loads,  resulting  from  the  de¬ 
formation  from  state  I  to  state  II  is  given  by  the  sum  of  the 
quantities  in  (2.15)  and  (2.21).  With  (2.11)  and  (2.12),  G[u] 
can  be  expressed  by 

G[u]  =  2^3  I  I  [(W^+a|w2){ic[u]  +  (A^+2A2)E^} 

+  a2w2{4A^(1-a2)Eu^^^  +  (a2-a2+3a![A^-3a2)e2 

+  2Ek[u]} 

+  ^{W^^i2+2A2w^2iCi*j)  *  ^3W22(i+j)^> 

+  W^^^  +  W^^^  -  2A^W^E^(1+E)“^ 

+  A^W2E^{i-E[A2(2+e)  -  2a2(1+e)-^ 

+  aJa^(4+3E)  (1+E)"^]}  +  W^^^JdCj^dCg  .  (2.22) 

We  shall  say  that  the  state  I  is  stable  if  G[u]  is  positive 
definite  for  all  u  lying  in  a  neighborhood  of  u  ■  0  and 
satisfying  the  kinematic  constraints  (2.5)  and  (2.19). 
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3.  Critical  equilibrivun  states 

A  necessary  condition  for  stability  of  state  I  is  that  the 
second  variation  Gg[u]  of  G[u]  be  non-negative  for  all 
sufficiently  small,  kinematically  admissible  values  of  u  .  With 
(2.10)  and  (2.12)  it  follows  from  (2.22)  that 


GgLu]  =  2^3  I  I  [k^{<[ul  +  (x2+2x2)e2} 

+  x2w2{4X^(1-X^)EU3^^^+A^e2} 

+  kgk^  +  2k23^A2Ek  +  , 

where 

ki  •  .  kj  -  2CWii*2x|k^2.X^W22)  - 

•'21-  > 

and 

k  =  X^(l-x2)u^^^+  (x|-x2)e  , 

Ai  =  (x2-x2)(3x2+^2j^-2  ^ 

Ag  =  -(x2-x2)(x2-x2)X-2  . 


(3.1) 


(3.2) 


(3.3) 


The  necessary  condition  for  stability  may  then  be  written  as 


GgEu]  >  0 


(3.4) 


for  all  u  satisfying  (2.19)  and  the  linearized  incompressi¬ 
bility  condition 


^1,1  * 


0 


(3.5) 
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which  is  obtained  from  (2.5). 

For  any  specified  value  of  ^3  »  ®  critical  value  of  A  , 

for  which  a  bifurcation  solution  of  the  static  problem  exists, 
occurs  when 


6G^[\i]  =  0  (3.6) 

for  some  non-trivial  displacement  field  u  which  satisfies  the 

boundary  conditions  of  the  problem.  We  shall  call  the  state  I 

corresponding  to  such  a  value  of  A  a  critical  state .  With  the 

definition  of  stability  given  at  the  end  of  §2,  state  I  will  be 

stable  if  G2[u]  is  positive  definite  in  a  neighborhood  of 

u  =  0  .  We  say  that  a  critical  state  is  at  the  stability  limit 

if  G„[u]  has  a  zero  stationary  value,  i.e.  if  G  [u]  =  0  for 

2  - 

a  non-trivial  value  of  u  satisfying  (3.6)  and  (3.5)  together 
with  the  boundary  conditions  of  the  problem. 

In  order  to  determine  such  a  value  of  u  ,  we  proceed  in  the 
following  manner.  We  take  account  of  the  constraint  condition 
(3.5)  by  introducing  the  Lagrange  multiplier  -4£^p(C^,C2) 
obtain  from  (3.1) 

SGgLu]  =  2^3  I  I  [k^{6K[u]  +  2(A^+2A^)E6E} 

+  a2w2{4A^(1-A^)(E6u^^^+u^^^6E)  +  2A^E6E} 

+  2k2k5k  +  2k23^A2(E6k+k6E) 

+  4A^W22a2e6E  - 

(3.7) 

where,  from  (2.10)  and  (3.3)^^, 


12. 


6<[u]  =  2{  [  (u^^^-Xu2^2^*^“i'^ 

1  ttU2_2-AUi_i)«U2* 

(“l,2*'“2,P%'' 

■  t“l,ll*“l,2S-'*“l‘ 

f^2,ll'*'“2.22^*“2^ 

k«k  =  x^Cl-X^)^{(Uj^_j^5Uj^) 

'l‘ 

+  X^(l-X^)  CX^-X^){(E6u^)  ,^+  "■ 


(3.8) 


With  (3.8)  and  (3.3)^^,  we  obtain  from  (3.7) 

f/2  -t, 

6G2[u]  =  4>^3|j  *^^2 

■^1  Z 

*  I  12* 


^  (3.9) 


where 


11 


21 


22 


33 


{k^.x2(l-x2)2k2)u^_^-Xk^U2_2 

♦  Xj^(l-X^)E{2x|w2*(x2-x2)k2*A2k2j^)-X‘’T>  . 

’  '’is  ’  ''l'“l,2*^“2.1>  • 

’'l<“2,2'^“l.l’  ■  ^2^’’  ’ 

B{k^a^*2\p  * 

*  CX?-X?)[Cx5-x|)k^*2A2k2j  *  SA^A^W^^) 


(3.10) 


■3  2^ 


+  X^(1-x2)U3^^^{2x|W2+  *  ^2^21^  ‘P  ' 

From  (2,19),  6u^  is  constant  on  each  o£  the  faces 

^1  “  *'^1  ■  '^hen,  from  (3.9)  the  condition  fiG^  =  0  implies 


"icc.  ■  "ac.a  ■  “  >  n  ''33«l'*«2  ■  “  - 

together  with  the  boundary  conditions 
/2 
-^2 

and 


^12  “  ^22  ■  ^  ^2  “  *  "^2 


Solutions  of  the  homogeneous  system  of  equations  (3.5) 
(3.11) -  (3.13)  can  be  found  of  the  form 


•sinflC- 


cosnCi 


U^(C2)  »  = 


cosQC^  ^ 


sinQC. 


>  U(?0  , 


E  »  0  ,  p  = 


cosfiC 


sinfiC- 


where 


n  =  nv/ZZ^  ,  n  =  1,2 , . . . , 


the  upper  (lower)  solution  corresponding  to  n  even  (odd) . 
The  functions  ^^(5^)  and  are  given  in  terms  of  U 


\Tl‘ 
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(3.11) 

(3.12) 

(3.13) 

and 

(3.14) 

(3.15) 

by 

(3.16) 


where  the  prime  denotes  differentiation  with  respect  to  Sg 
and 
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e  =  u"-{i+(x-i)^A}n^u  ,  (3.17) 

with 

CX  +  1)% 

^  '  ^^3^1  '  ’  (3.18) 

U(C2)  satisfies  the  differential  equation 

u^^'^^-{x^+i+(x-i)^A}n^u"  +  x^nS  =  o  ,  (3.19) 


and  the  boundary  conditions 


u"  +  x^n^u  =  0 

u'” -{2X^  +  l+(X-l)^A}fi^U'  =  0 


(3.20) 


Solutions  of  (3.19)  with  (3.20)  have  already  been  obtained 
(see,  for  example,  [1,2]).  These  are  either  even  or  odd  func¬ 
tions  of  ^2  represent  respectively  flexural  or  barreling 

deformations . 

From  (3.1)  we  obtain,  with  (2.10),  (3.3)^,  (3.10)  and  (3.14)2, 


G2[u]  =  2^2  I  I  {k^K[u]  +  x2(l-x2)2k2u2^^}dCid^2 


(3.21) 


We  have,  with  (3.12)  and  (3.13), 

1.  o  1. 


[  =  j  C"aPa2]  '  |  [u^F^J  "^1^2  =  0 

~^2  "n 


(3.22) 


15. 


Then,  with  (3.11)^^  2»  (3*5),  (3.14)j  and  (3.22),  equation  (3.21) 
yields 

GgLu]  =  0  .  (3.23) 

Accordingly,  for  any  displacement  field  of  the  form  (3.14)  the 
critical  states  are  states  of  neutral  stability. 

Substituting  from  (3.14)  in  (2.10)2»  we  obtain,  with 
(3.16)^, 

k[u]  =  (2A^(2^)"^[(A^+l)fl^U'^+  U"^+  A^£2‘*U^+  2A^jJ^(UU')' 

+  (-l)'^cos  2flCi((X^+l)G^U’^-u"^-A^n^U^-2A^n^a}]  ,(3.24) 

where 

a  =  UU"  -  U'^  .  (3.25) 

We  note  that  k[u]  is  necessarily  an  even  function  of  Cg  . 


L 


X 


4 .  Potential  energy  near  a  critical  state 

Paralleling  the  procedure  of  our  previous  paper  [3],  we 
develop  the  stability  condition  that  G[u]  be  positive  definite 
for  all  u  satisfying  the  kinematic  constraints  (2.5)  and  (2.19) 
and  lying  in  a  neighborhood  of  u  *  0  . 

We  take 

u  *  eu  +  e^u  ,  (4.1) 

where  u  is  the  solution  for  u  given  by  equations  (3 . 14) - (3 . 20) , 
e  is  a  small  parameter,  and  u  satisfies  the  orthogonality  con¬ 
dition 


We  then  determine  the  value  of  u  which,  for  a  fixed  value  of 
e  ,  gives  a  stationary  value  of  G[u]  .  If  this  value  of 
G[u]  is  positive  (negative)  then  the  critical  state  is  stable 
(unstable) . 

With  (3.14) j  and  (2.4),  we  can  rewrite  (4.1)  as 

“d  *  ^ 


where,  from  (4.2), 


0  . 


(4.4) 


Also,  since  u  must  satisfy  the  incompressibility  constraint 
(2.5),  we  obtain,  with  (4.3)  and  (3.5)  and  the  neglect  of  terms 


17. 


of  higher  degree  than  the  first  in  e  , 


2  “  ^  2^2  1~^1  1^2  2^ 


Again,  since  u  must  satisfy  the  constraint  (2.19)  on 
5^  *  ±  ,  we  have 


where  e  is  a  constant. 

It  is  shown  in  Appendix  A  that,  with  (4.3),  G[u],  given  by 

(2.22),  may  be  expressed  to  order  in  the  form  (A8) ,  thus 


G[u]  =  G[eu+£^u] 

=  e^G2[u]  +  2l^e^  j  |  (gp^+gp^+eg^^^dS^d^g  ' 


where  the  g's  are  defined  in  (A9) .  From  (3.23),  G2[u]  =  0. 

We  will  now  show  that  the  remaining  terms  on  the  right-hand  side 
of  (4.7)  are  of  order  . 

We  replace  Ujj,p,E  by  u^,p,E  »  0  in  (3.10)  and  use  the 
resulting  equations,  together  with  (A3)  and  (A9)j^,  to  express 
gp^  in  the  form 


(4.8) 


From  (3.14)j^  2»  (3.16)j^and  (3.25)  we  obtain 


18. 


“l,2“2,l  -  “l.A.S  ■  (-l)"ccos2a5^)  .  (4.9) 


Also,  from  (3.14)^^  (3.15),  (3.24)  and  (4.9),  we  obtain 

ll 

I  ■  I  •  I 


-t 


-t 


t 

* 

-I 


We  now  substitute  from  (4.5)  in  (4.8).  From  the  resulting 
expression  for  and  (4.10),  (3.11) -  (3.13) ,  we  obtain,  to 

order  e  , 

I  I  -  2eX3  [  f 


A1 


so,  from  the  expression  (A9)2  for  and  (4.10),  we  obtain 


I  I  gp^dC,d5^  -  0  . 


(4.12) 


With  (4.11),  (4.12)  and  ^2^^!  “  ®  »  equation  (4.7)  yields 


G[u]  =  e^G[u]  =  2^3e'♦  j  |  gdC^dCg  . 


where 


(4.13) 


g  =  g 


(4.14) 


T 


5 .  Conditions  for  G  to  have  a  stationary  value 

In  this  section  we  shall  obtain  differential  equations  and 
boundary  conditions  determining  the  field  u  which  satisfies  the 
kinematic  constraint  (4.5)  and  the  orthogonality  constraint  (4.4) 
and  leads  to  a  stationary  value  for  G  .  We  use  the  method  of  un¬ 
determined  multipliers  to  remove  the  constraints  and  accordingly 
write,  neglecting  terms  of  order  e  in  (4.5), 


G[u]  .  2i3  J  j  (g  -  ® 

^  ^>‘“a,6“a,8J‘‘5l'‘«2  "  f^.l) 

where  and  At^x  =  constant  are  the  Lagrange 

multipliers  associated  with  the  constraints  (4.5)  and  (4.4)  res¬ 
pectively.  Then  for  G[u]  to  have  a  stationary  value. 


6G[u]  =  0 


(5.2) 


for  all  u  in  the  neighborhood  of  u  =  0  . 

From  (5.1)  we  obtain,  with  (4.14),  (A9)4,  (Al)  and  (A3), 


6G[u]  = 


«3  I 

-  I  1  ffa8,6S-?'33®««l««2}' 


(5.3) 


where  F  ^  and  E,-  are  defined  by 

Otp  JJ 


20. 


-X^CX^-l){2x2w2-(x|-x2)k2+A2k2i}E 
-XiCX^-l)k2{^[ul  ^  U^,i(Ui^1-Au2^2^^ 


-12x3(x2-l)^k, 

if“i.i>  *  ’‘"1,1 

> 

^21  = 

*  HP“i,2 

-X^(X^-l)k2U^^ 

><“2,1  - 

^12 

*  HP"2,1 

-X^(x2-l)k2Ui^ 

><“1,2  ’ 

II 

OJ 

CM 

Itl4 

^l(“2.2-^^l.l) 

-  x‘^p-XjJS^_^ 

-X3^(X^-l)k2U^ 

.i^“a.2‘^“i,i^  * 

X1^2,2  * 

F33  =  {(x|+2x|)k^+x|Aj^W2+(x2-x2)2k2 
-2(X=-xf)A,k2,*2X^A2^,,}E  -  p 
-X,(x2-l){2x|l»3-(x2-x|)k^*A^lc^^)u^_^ 
lCX^x2)k^. 

+  2x2(x2.1)2{k2^.6(x2-Xpk3+2A2k3^}(u^^^)2. 
With  (5.3)  and  (4.6)  the  condition  (5.2)  yields 

^la,a  *  ^2a,a  *  |  f  ° 

S  - 

^21  *  I  ^1^^2  “  0  ^1'  "^1  • 


and 


(5.4) 


(5.5) 


(5.6) 


F  »  F  »  0 
12  22 


on  ^2“  *^2  • 


21. 


These  relations  have  the  same  form  as  (3 .11) -  (3. 13)  with 
replaced  by  F  «  . 


Neglecting  terms  o£  order  e  in  the  constraint  (4.5),  we 


obtain 


+  X"J’‘  u 


2  2,2 


-E  +  X-(u  U-  -u  u  ) 

j  Jl  f  ^  X^X  ^  ^ 


(5.7) 


With  the  expressions  (5.4)  for  F^^g  and  F^^  and  the  expressions 
for  u  and  p  obtained  in  S3,  equations  (5. 5) -(5. 7),  together 
with  the  orthogonality  condition  (4.4)  and  the  additional  boundary 
condition  (4.6),  provide  a  set  of  equations  for  the  determination 
of  u^,  E,  e,  p  and  x  • 


6. 


Development  of  the  governing  equations  for  u 
In  this  section  we  obtain  from  (5. 4) -(5. 7)  expressions  for 
u  analogous  to  those  given  for  u  in  (3.14). 

We  first  use  the  expressions  obtained  in  §3  to  substitute 

^  Ak  —  ^ 

for  u^  and  p  in  the  expressions  (5.4)  for  and  F33 

and  obtain  the  new  expressions  for  F^^  and  F^^  given  by 
equations  (B1)-CB3)  of  Appendix  B. 

From  (3.14),  (3.15)  and  (4.6),  we  obtain 


u 


1,2 


u 


2,1 


u 


1,2 


on 


£  =  +  £ 
^1  -  ^1 


(6.1) 


Accordingly,  with  the  expressions  (5.4)2  ^21’  boundary 

condition  (5.6)^^  may  be  replaced  by 


Ug,!  =  0  on  5^*  ± .  (6.2) 

From  (5.7)  and  (4.9)  we  obtain 

=  (2X2)"^{ (UU')'-(-l)^acos2n£^}  .  (6.3) 

We  shall  assume  a  solution  for  u^ ,  E,  e,  p  and  x 
equations  (5. 5) -(5.7),  (4.4)  and  (4.6)  of  the  form 


u^  »  U^sin2n£^  +  » 

u  =  U  cos2D£  +  ^  UU’-X  (e+E)£„  , 
^  1  2  2 

p  *  P  cos2n£^  ♦  Q  , 

X  =  0  , 


where  U,  U^,  P  and  Q  are  functions  of  Cg  o’^ly,  U  is  a 


23. 


function  of  Cg  only  determined  by  equations  (3.19)  and  (3.20). 

We  note  from  (3.14)  and  (3.15)  that  expressions  for  u^ 
of  the  form  (6.4)^^  2  automatically  satisfy  the  orthogonality  con¬ 
dition  (4.4),  the  boundary  condition  (4.6),  and  the  boundary  con¬ 
dition  (6.2)  which  is  equivalent  to  (5.6)j^.  By  substituting  from 

(6.4) i  2  ill  (6*3)  it  follows  that 

where  a  is  defined  in  (3.25). 

Introducing  the  expressions  (Bl)  for  given  in  Appendix 

B  into  (5.5)^  2»  obtain  with  (6.4)^  and  (3.18)  the  differential 
equations 

.  c-i)”ui2-mn)sin2ae^  , 

where  fg^  and  (|)^g  are  defined  in  equations  (B2)  and  (B3) 
of  Appendix  B.  Now,  substituting  in  (6.6)  for  u^j  and  p  from 

(6.4)  12  3  obtain 

k^{ir^-4(22[l+(A-l)2A]U^}  +  2nA-£^P  »  (-l)’'((f^2-2««3^^)  , 
k^(u”-4n2u)  -X‘lp’  -  (-l)"(4»22+2n(fr23^)  ,  (6.7) 

k^(UU’)”  -  2Q’  -  2X2f’2  . 

Eliminating  P  from  (6.7)^.  2  using  (6.5),  we  obtain 

{X^+1  +(X-l)^A}ir'  +  16X^Ji**U  -  -(-1)'^X^4'  ,  (6.8) 


24. 


where 


<&’  =  X-^{q'  *  ^  C4,22+2n(t>2^)}  . 


and 


1  n  o  p  4AXqJ2 


(6.9) 


It  is  shown  in  Appendix  B  that  $  may  be  expressed  in  the  form 
given  in  equation  (Bll). 

With  (6.5),  equations  (6.7)^^  2  yield 

—  X  k  _  o  — 

P  =  {  U'"  -4fl^[l+(X-l)^A]u’  +  (-l)^q}  , 

4Xn‘=' 


(6.10) 


Q 


K(UU’) 


2  22 


With  equations  (Bl)  in  Appendix  B  and  (6.4)^,  the  boundary 
conditions  (5.6)^  ^  yield 

^1  ^^1,2'^ =  (-I)%^2sin2n5i  , 

on  ^2"  -^2  •  (6.11) 

Then,  substituting  from  (6.4)  in  (6.11),  we  obtain,  with 
(6.5)  and  (6.10)^^  , 

U”  +  =  -(-1)^X^$^  , 

on  Cg  =  ±  ^2  (6.12) 

iJ'” -4J2^{2X^+  CX-1)^A}U'=  (-1)'^X^$2 


25. 


where 


= 


= 


,  4xx,n 

2X„ 


2x7^  -1 


~T^  *22^^ 


(6.13) 


It  is  shown  in  Appendix  B  that  and  9^  may  be  expressed 

in  the  forms  given  in  equations  (B14). 

We  also  obtain  on  ^2  "  *  '^2* 

Q  =  k^[|cuu')'-x^(2e  +  E)]-X2£22  •  (6.14) 


It  follows,  with  (6.10)2.  that  (6.14)  is  valid  throughout  the 
body. 

The  relations  (5.6)2  (5.5)^  express  the  assumptions  of 

dead-loading  in  the  1  and  3  directions  respectively.  We  sub¬ 
stitute  in  them  the  expressions  for  and  given  in 

equations  (Bl)  of  Appendix  B.  Then,  using  (3.2),  (3.3),  (3.15), 
(3.18),  (6.4),  (6.5),  (6.9),  (6.10),  and  (6.14)2>  we  obtain 


'^e  +  aE  *  X^b^  ,  ae  +  =  X^b^  , 


(6.15) 


where 


a,  = 


3X^  +  1  +  (X-l)^A  , 


a 


*1*1 


(6.16) 


26. 


1 


(-D^'X  rt  , 

- 5-i 

U.Sl^  J 

'P 

+  ^  f  2  {  mn '  •)  '  +  -2- 


^  {(UU’)'  +  -~  [1+(X-1)2a]  + 
-  2>2 


ik-J  '  <C™')'  *  ^  tf33-Va2»'*52  • 


(6.17) 


With  the  relations  (6.12)^^,  (6.13)^  and  (6.9)2»  equation 
(6.17)^  yields 

\  ^  {(UU’)'  +  (f^^-Xf22)}d52  •  = 

Co  X 


i 6 . 18) 


The  expressions  for  and  may  be  rewritten  in  the  forms 

given  in  equations  (B16)  of  Appendix  B. 


Equations  (6.15)  may  be  solved  for  e  and  to  yield 
-  ^3  -  ^3 

®  *  X  Cb^a^-b^a)  ,  E  =  x  ’ 


(6.19) 


where  h  ,  assumed  to  be  non- zero,  is  defined  by 


A  = 


From  (6.19),  we  obtain 


biC  +  bgE 


(a,b?+a  bj-2ab  b^) 
A  21  12  12 


(6.20) 


(^.21) 


27. 


7 .  Development  of  the  expression  for  G[u] 

In  this  section  the  expressions  obtained  for  u  and  p  in 

m0 

§3  and  the  expressions  obtained  for  u  in  §6  are  used  to  write 
the  expression  for  G[u]  given  by  (4.13)  and  (4.14)  in  terms  of 
the  functions  U(C2)  and  UCC^)  only.  With  (2.10)2*  (5.4), 
(5.6)j^  and  equations  (Al) ,  (A3)  and  (AO)^  in  Appendix  A,  we  can 
rewrite  the  expression  (4.14)  for  g  as 

g 

where 


«11  ■ 

-AjPU 

2,2' 

X^(x2.1)k2{|<[u]+u^ 

-12xJ(x2 

A  A 

-X3PU 

l.l' 

X^(X  ■1)^2^1,1^“2 ,2 

-XUi,i)  , 

Hi2 

A  A 

X3PU2 

,1-^ 

l(X^'l)k2“i,i(^l,2'^ 

XU2^l)  , 

»21  = 

A  A 

X3PP1 

,2-' 

»33  - 

<»3*2 

4(X 

=-x2)lt24A2k2i)K[ul 

H  =  ik2C<[u])2+l2x2(x2.i)2],^((;^^^)2^[^] 

+  16xJ(x2.l)‘»k^(G^^^)‘'  . 

We  now  substitute  in  (7.2)  the  expressions  for  u^  and  p 
given  in  (3.14)  and  use  (3.16),  (3.18)  and  (3.24).  Then,  with 
the  notation  introduced  in  equations  (B4)2»  (B5)  and  (B17)  of 


+  X3p(u^  2^2  1^2  2^  * 


(7.1) 


*•'  T.  «*•.. 


28. 


Appendix  B  and  the  further  notation  Cc£.(3.24)) 

+  u"^  +  +  2A^fi^(UU')’  , 

<2  “  +  -  ZX^Sl^a  , 


(7.3) 


h^l  =  2b’u'  +  A{ic^+2(x2+l)n^U'2}  +  4(X-l)^Bn2u’2  , 

hgg  =  b'u'  +  2X^AJl^u’^  ,  (7*5) 

h  -  CK^  +  (x^-i)^  dn^u'^  , 

33 

and 

-  2e'u'  +  A{iC2+2(X^+l)n^u'^}  +  4(X-l)^Bn2u’2  , 

’^22  *  ^22  ’ 

^33  "  ^'^2  *  » 


(7.6) 


29. 


ni2  =  3'u  +  Au' (u"+A^fl^U)  , 

Hgi  =  p'u"  +  x^An^u* (u"+n^u)  . 


We  also  obtain 

k 

H  {XX2h+C-l)^h^cos2nC3^+h2Cos4n5^}  ,  (7.7) 

where 

h  =  - 1 -  (2ic^+»c|)  +  n^u'^(2K  +K  ) 

4CX^-1)  ^  ^  ^  ^ 

+  (X-l)'*Cn‘^u''*  ,  (7.8) 

and 

48(X+1)H, 

n  —  _ _ _ _ _ _  *♦ 


The  corresponding  expressions  for  h^  and  h^  will  not  be 
required. 

We  now  substitute  from  (7.1)  in  (4.13)  and  use  (4.4),  (5.5)  and 
(5.6)  to  obtain 


G[u]  =  2f, 


f  \ 


B  *  ®«33  *  » 


(7.9) 


Then,  introducing  into  (7.9)  the  expressions  (7.4)  and  (7.7)  for 
and  H  and  the  expressions  for  u^  and  p  given  by  (6.4) 
and  (6.5),  using  (4.9),  and  carrying  out  the  integration  with 
respect  to  , 


we  obtain 


30. 


G[u] 


{(-l)"X2g^+g2}dC2 


Ceg3+Egj^)dC2+ 


where 


(7.10) 


=  -aiT”  +  2113^2^*  {4n^a+4(X^-l)Af2^a  +  2t)^^ 

*  4n22>u'  +  Sn^^U  , 

82  '  ’^11“  ^  ^12“'  ^hggCUU’)’  -X2aq  ^ 

.  4xV{(UU')’>2  -  CUU')',  (7.11) 

>=1 

+  2h22  +  4XV(UU')’  , 
gj^  »  2{h22  +  X^n^fuu’)’-  h33)}  . 


Introducing  into  (7.11),  .  the  expressions  (7.5)  for 

0,4 

^11’  ^22’  ^33»  obtain 

g3  =  46'u'  +  4X^n^(UU')'  +  2(3X^+l)An^U’2  +  Ak^ 

+  4(X-l)^Bn2u’^  ,  (7.12) 

gj^  =  2{b'u’+  X^fl^(UU')'-  c<^+  2x^An^u'^ 

-(X^-l)2dfi^u'^}  . 

From  the  expressions  (B7)3and  (7.3)  for  g*  and  we 

obtain,  with  (B15)^  , 


31. 


3'u’  =  (U'u")'-u"2  -{l+(X^-l)A}n^u'2  , 

=  [u' cu"  +  x^n^u)  -  U{U"'  -[2X^+l+CX^-l)A]fi^u'}]  '  (7.13) 

-(X^-l)AJ2^u'^  . 


We  now  substitute  from  (7.13)  in  (7.12),  use  (BIS),  and 
integrate  the  resulting  expressions  for  g^  and  gj^  .  Then  with 
relations  (3.20)  we  obtain 

i  |^^(Fg3+Egi^)d?2=  2£2Afl^(b^e+b2E)  ,  (7.14) 

-^2 

where  b^  and  are  given  by  (B16) . 

Noting  from  (7.6)^,  with  (B7)3,  that  n2i  be  written 

in  the  form 


n 


21 


1 

7 


n  2i» 


where 


^21  “  U"2  +  (A-l)fi2u'2+Ax2n**u2  , 

and  that 


air"  -  (air*) '-(a'u')  '  ♦  a"u'  , 


_  r  77*  ^  *  '  77' 

n^2U  -  (n^2U  )  -n  , 


(7.15) 


(7.16) 


(7.17) 


we  see  that  the  expression  (7.11)^^  for  g^  may  be  rewritten 
in  the  form 


32. 


®1  “  *  ®12» 

where 

®12  “  ■“^'  (a'+2n^2^^'  ‘ 


(7.18) 


(7.19) 


The  expression  (7.19)2  may  be  rewritten  (cf.(6.12)2)  in  the  form 

gl2  =  -a{U'’+4X^f2^U+(-l)^X^*^} 

+  (a'+2ni2)U‘+4(X^n^o+W2i)U+(-l)”x^a*3^  .  (7.20) 


Using  (6.12)^,  (6.13)^,  (7.6)4,  (7-16),  (Bll) ,  (B13)^  and 
(3.20),  we  obtain  from  (7.20) 


\  ^  ®12‘^^2  “  [{7x2+i+(x2-i)A}n2uu'u' 


-4n^{ (X^+l-A)U'^-x2An^U^}U]  ^ 

~^2 

n^{7x2+i+(x^-i)A}[uu'  (u'^+x^n^u^)]^^ 

2 


(7.21) 


In  order  to  simplify  the  expression  for  "f®  substi¬ 

tute  from  (7.6)  and  (7.16)  in  (7.19)^^  and  employ  (B7)  and 
(7.3)2  to  obtain 

-  -36(u"  +  X^n^U)  +  68'u'-A{6X^J2^UU"+4u"^ 

+  2u'u'"+  3(X^+l)X^n^U^-6(2x2  +  l)n^U'^} 

+  8(X-l)^Bfi^U'^  , 


(7.22) 


33. 


where  B  is  given  by  (B7)^. 

From  (7.11)2  we  obtain,  after  a  lengthy  calculation  in  which 
equations  (82)2,  CB3)^^5,  (84)2,  CB.5),  (6.9),  (6.10)^,  (6.14), 
(7.3)2t  (7.5)2,  and  (7.6)^  are  employed, 

%2  ~  ’  (7.23) 

where 


=  a(a’-26’u)  +  2AaU ' (u"+X^n^U)  , 

and  (7.24) 

822  “  4n^a^{l+(X^-l)A}  +  4a(6’u’-6"U) 

+  166'U'(UU')'  +  8X^n^{(UU’) ’}^ 

+  4A{a(3(x2  +  l)n2u’^-U’'2-x2n’*u2-2x2a2al 
+  a’u’ (U"+X^Q^U)  +  8x2n^U’2(UU’) '} 

+  16(X-l)38n^aU'2  . 


With  (3.20)j^,  (7.24)^  and  (812),  equation  (7.23)  yields 


-I 


■^^(x^-i) (1-A) [uu’ (u’^+x^n^u^)]  ^ 

~^2 


T  I  ^22^^2  • 


(7.25) 


We  now  substitute  in  (7.10)  from  (7.11),  (7.18),  (7.21), 
(7.25)  and  (6.21)  to  obtain 


G[u]  -  4^^£2£3X3k^(G^-G2+G3-G^)  , 


(7.26) 


34. 


where 


G,  = 


1  f^2  r  2X^  „  ,7.^4 

j  2»  J  [§22^^”^^  SttG  ^  h]d^^  , 


2  -Z. 


X^ 
^3 


S  [uu'{u-2.xW}]^^^,^  , 


■^'/^  [{7x2  +  1+(X^-1)A}UU'u' 

8  X  X^Z2 


G.  = 


(7. 


-4u{u'^(x^+i-s:)-jx^n^u^}]  „  , 

^2"  2 

A  (aib2.a2b2-2abib2)  . 


2^2„2- 


G,.  = 


where  A  is  defined  in  (6.20)  and  the  a's  and  b's  in  (6.16) 
(6.17)2  (6.18), 


27) 


In  this  and  the  following  sections,  we  suppose  that  the 
plate  is  thin  and  calculate  G[u]  for  the  case  when  the  bifur¬ 
cations  is  of  the  flexural  type. 

It  has  been  shown  in  [  3  ]  that  if  the  value  of  n  de¬ 
fined  by  (cf.C3.15)) 


n  =  0^2  " 


(8.1) 


is  small,  then  the  critical  value  of  X  at  which  a  flexural 
bifurcation*  can  occur  is  given  by 


X-1  =  +  O(n^) 


(8.2) 


We  now  introduce  the  notation 


^(0)  ^ 

’  a 

X«=l 

w ( O)  _  M  1 

’  aB  ■  o3 

II 

r< 

etc. , 

^a|x»l 

k(o) 

’  o8 

^obIx*! 

t,8=l, 

,2)  , 

(2-9)i,2 

and  (2. 

10)j^  that 

*  x-i(x 

-1)2{1-CX-1) 

+  (X- 

■1)^ 

-(X-l)^+  0(X- 

•dS 

(8.3) 


(8.4) 


•  The  integer  n  in  (8.1)  is  the  number  of  half-wavelengths  along 
the  1-direction  in  the  flexural  mode  under  consideration.  If  the  lowest- 
order  mode  is  not  suppressed  hy  means  of  sane  type  of  passive  constraint,  no 
value  of  n  greater  than  unity  is  relevant. 
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With  (8.3)  and  (8.4),  Taylor’s  theorem  yields  the  following 
approximate  expressions  for  W^,  W^^,  ^aey6  ' 

=  W^°^  +  X-l(X-l)2[l-(X-l)  +  (X-l)2](wf°Ux2w^°^) 

+  2^3  (A-lJ  C\ia*^^3^12cx*^3'^22J 
+  0(X-1)5  ,  (8.5) 

"aey-  ’ 

”aers  •  "aer-s*  ' 


From  (8.3),  (3.2)^^  2  (A7)j^  we  have 


W^Q^+X^W^^^  ,  k^°)=2(wJj^+2X^wf°)+X^wf°)) 

*=3°’  '  • 


(8.6) 


With  (8.5)  and  (8.6),  equations  (3.2)^^  2  (A7)j^  ^  yield 

ki  =  k^°^+jX"^k^^^(X-l)^[l-(X-l)  +  (X-l)2]+3X“^k^°^(X-l)'^ 

+  0(X-1)5  , 

(8.7) 

kg  »  k^°^  +  12X2^kf°^(X-l)2[i-(X-l)]+0(X-l)'*  , 

k3  =  k^0^+4X3^k^°^(X-l)^+O(X-l)3,  k^»kj°^  +  0(X-1)2  . 


From  (7.8),  (B4)2,  and  (B5) ,  we  obtain  with  (8.7) 


A  -  +  (2-A^°^)  +  B^°^}(X-1)^[1-(X-1)]  +  0(X-1)‘*  , 

(8.8) 

B  -  B^°^-yB^°^ (X-1)  +  O(X-l)^,  C  • 


+  O(X-l)  , 


37. 


where,  from  (B4)2,  (BS)  and  (A7)^  , 


=  A| 

X=1 

4kf°V(X3k{°b  , 

=  B| 

II 

H 

II 

48k^°V(X3kp^), 

C8.9) 

C  =C  1  = 

32 

/•ui(O)  +4X^W^^^ 

■  ^”1111  ^^3^1112  °^3i122  ^^3^1222 

^  ^|x=l 

A3W2222J 


I 

I 


i 

t 

! 
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9 .  Asymptotic  first-order  solution 

We  define  the  dimensionless  thickness  coordinate  t  by 


t  = 


C9.1) 


and  write 


V(t)  *  0(^2)  . 


Then,  noting  from  (3.15)  and  (8.1)  that 


(9.2) 


Id  =  1  d 

n  352  ^  ^ 


and  using  the  expression  (8.8)^  for  ,  we  can  rewrice  the 

differential  equation  (3.19)  as 


(9.3) 


-4^  -{2+2(X-l)  +  (X-1)2(1+A^°4  +  0(X-l)^}-4 
n  dt 


+  {1+2(X-1)  +  (X-1)^}V  =  0  , 


(9.4) 


and  the  boundary  conditions  (3.20)  as 


2 

1  IV  +  {1  +  2(X-1)  +  (X-1)2}V  =  0  , 
dt^ 


^  ^  -{3+4(X-1)  +  (X-1)2(2+A^°4  +  0(X-1)S  ^  ^  =  0  , 

r,3  ^^3  n  dt 


(9.5) 


when  t  =  ±  1  . 

With  the  asymptotic  expression  for  X-1  given  in  (8.2), 
equation  (9.4)  can  be  rewritten  as 
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+  {1 +|n^  +||n‘*+0(n^)}V  =  0 


(9.6) 


and  the  boundary  conditions  (9.5)  can  be  rewritten  as 


40^  {1.4,2,  52^4  ^ 

n  at 


1  ^  -  {3.|n2.4,4(^,^(0)3  ,  0(n^)}^^=  o 

dt^  3  y  b  dt 


(9.7) 


on  t  =  ±1  .  It  can  easily  be  verified  that  the  solution  of 
(9.6),  with  the  boundary  conditions  (9.7),  is  given  by* 


V  =  1- 


jH  t  +n  (jt  -^t  )  +  n  (jj-t  -jg-t  )+  o(n  )  .  (9.8) 


We  note  that  to  order  ,  V  is  independent  of  both 
and  the  form  of  the  strain-energy  function  W  . 


*  We  employ  the  normalizing  condition  V(0)  =  1. 


40. 


(10.4) 

and  tCCg)  »  given  by  (Bll)  and 

(6.13) . 

With  (9.1),  (8.2),  (8.8)  and  (B5)  we  obtain  from  (10.4), 


(Bll)  and  (B14) 
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^  =  4nn2t(2+A*^°b  +  O(n^)  , 

-*i|t=±l  "  ±2nn[2+^n^  +  n‘'(^  +  +  0(n^)]  ,  (10.5) 

^2|t=±i  “  2nn^[i+n^(2+^A(°b  +  oCnS]  . 


With  these  expressions  introduced  on  the  right-hand  sides  of 
(10.2)  and  (10.3),  it  can  easily  be  verified  that  the  solution 
of  (10.2),  subject  to  the  boundary  conditions  (10.3),  is 


^  +  +0(n9)  , 

(10.6) 

where  D  is  a  constant  which  could  be  evaluated  if 
the  calculations  were  carried  out  to  a  higher  order  in  n  . 
However,  the  value  of  D  will  not  be  required  for  the  calculation 
of  G[u]  which  is  the  main  object  of  the  present  paper. 

It  is  seen  from  (10.6)  that  V  -  and  hence,  from  (10.1), 

U  -  depends  on  the  form  of  the  strain-energy  function  W  only 
through  A^®^  ,  defined  in  (8.9)j^.  From  (8.6)j^  and  (8.9)j^,  we 
see  that 


a(o) 


8(W^J^+2A2w^°^+A^wf°^; 

m 


m: 


X3(Wi 


X^Wg 


') 


(10.7) 
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11.  The  asymptotic  expression  for  G[u] 

In  order  to  calculate  G[u]  in  the  asymptotic  case  when 
1  »  we  use  the  formulae  (7.26)  and  (7.27)  and  introduce 
into  the  latter  the  relations  (9.1),  (9.2),  (9.3)  and  (10.1). 
Then,  we  introduce> in  the  expressions  for  G^,  G^,  G^  and  G^ 
so  obtained> the  expressions  (9.8)  and  (10.6)  for  V  and  V  and 
systematically  neglect  terms  of  higher  degree  than  the  fourth  in 
n  . 

Calculation  of  G^ 


Following  this  procedure  we  obtain  from  (7.27)^^, 


p  _  Q 
- 


32X- 


I  ^^2*^  Kl  ' 


(11.1) 


-1 


where  g^^,  g^^  and  h  are  defined  by 


g 


11 


11 


^  6 
§22  ""  ^  ^22  ^ 


h  =  n*®h 


(11.2) 


and  gii,  g22  ^  defined  by  equations  (7.22),  (7.24)2 

and  (7.8)j^.  With  (9.2),  (9.3)  and  (10.1)  and  the  notation 

a  =  Q-^a  ,  g  =  ,  <2  = 


we  obtain  from  these  equations 


;  .  ,  1  d-^V  .  .  6  de  dV 

>11  ■  36(-j  t-2  *  ^  *  “2  at  at 

n  dt  n 

jf  6X^  w  d^V  .  4  .  2  dV  d^V 


43. 


3(x2+1)xV  -A(2x2+l)(g:)2j 


+  -L(X-1)3b(^)2  , 
n 


’22 


,dV  d6 

„_(3t  3t-V’ 


'>7 

^  16  de  dv  d  ^  SX^fd  rvdv;\2 

*  nr  3t  3t  ^  nrlst 

.  4a(S[4cx2.1)(^)2  .  Jj_(d?V)2  .  xV  .  2x2^] 
1  da  dV,d^V  .  ,2^2v-i  .  8X^,dV<,2  d 

*  nr  It  ^  It^^iT^j 


*  ^(x-D^Bacg)^ 


h  = 


4(X‘'-1) 


*  !a^  = 


Tl"(X+l) 


*  XA4I 

4 

T1 


ccg)** 


From  (11.3)  and  the  expressions  (3.25),  (3.7),  (B.5)  and  (7.3 
obtain,  with  (9.2),  (9.3)  and  (10.1), 


a 


•jei-*') 

n"  dt2  I 


x2.i 


l+(x2-l)Ajv  , 


,dV^2  ^  1  fd^V.2^  2„2  ^  2X^  d  ,„dV^ 

^it^  *  nr^T^J  :X'  at  ^  w  ’ 

n  dt  n 


X‘^+l,dV.2 


1  ,d‘'V^2  .  .2„2  ,,2;; 

XV  -  2X  a 


From  (B5),  (8.8)j^  and  (8.2),  we  obtain  the  following 
asymptotic  expression  for  A  : 


(11.4) 


,  we 


(11.5) 
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A  =  +  O(n^)  .  (11.6) 

We  now  substitute  in  (11.5)  the  asymptotic  expressions  for 
V,  X  and  A  ,  given  in  (9.8),  (8.2)  and  (11.6)  respectively, 
and  obtain  the  following  asymptotic  expressions  for  a,  3,  ic^ 
and  <2  • 


a  =  -i+n^(|-2t2)  +Ti*‘(^-|tS  +  o(n^)  , 
g  =  -2+n2(|- t^)  +  r,'*(^-^Af°)  - -^t^)  ♦  o(n^)  , 

“  n^(-^+4t2)  +  n^(^- st^  +  8t'^)  +  o(n^)  , 

With  (9.8),  (8.2),  (11.6)  and  (11.7),  equations  (11.4)  yield 


=  n^(12-^A^°b  +  n‘^{4- ll-A^®^  +  4(2+A^°bt^}  +  O(n^)  , 

^22  "  4-Ti2(8+24t2)  -  n'*{||. -^(l-A^°))t^+40t^} 

+  O(n^)  , 

h  =  (j- jt^  +  3t**)  +  O(n^)  . 


Substituting  from  (11.8)  and  (8.2)  in  (11.1)  and  carrying 
out  the  integration,  we  obtain  the  asymptotic  expression  for  G^: 


.n2(|+^A^°b  +  +  0(n^)}.  (ll. 


TM  360 


.9) 


Calculation  of  G^ 

With  (9.2),  (9.3)  and  (8.1)  the  expression  (7.27)2  ^2 


45. 


rewritten  as 


J  't»i  • 


(11.10) 


With  the  asymptotic  expressions  (8.2)  and  (9.8)  for  X  and  V 
respectively,  this  yields 


C  —  16U-.  ^  ^  \ 


(11.11) 


Calculation  of 

With  (9.2),  (9.3)  and  (10.1)  the  expression  (7.27)2  for 
can  be  rewritten  as 


■  -57V 


8X- 


-4{(x2  +  1-A)  (^)^  >  Ax2nV}V]^^^ 


(11.12) 


We  introduce  the  asymptotic  expressions  (9.8),  (10.6)  and  (8.2) 
for  V,  V  and  X  respectively  and  obtain,  with  (11.6), 


n^(- 


U  l4^  •' 


n'*f25  - 


3F0^  ^ 


0(T1^)} 


(11.13) 


Calculation  of  Gj^ 

With  (9.1),  (9.2)  and  (9.3),  we  can  rewrite  the  expressions 


(B16)  for  b^  and  bg  as 


46. 


‘’l  = 


^2  = 


’IT 


n  at  n 

+  ^(X^-DA^  -  CX-l)^B}]dt  , 
^  n  at  n 


C11.14) 


-  [X^+l+CX^-l)c]A+CX''-l)''a‘}]dt  , 


where  c  and  d  are  defined  in  (B.17).  We  introduce  the 
asymptotic  expressions  for  X,  V  and  A  given  in  (8.2),  (9.8) 
and  (11.6)  and  obtain 

=  2b2  =  -fl^{l-|n^  +  +  0(n^)}  .  (11.15) 

Substituting  from  (11.15)  in  (7.27)^  and  using  (6.20),  we  obtain 

Gi^  -  2'FSi^{l-2Ti^ +  +  0(n^)}  ,  (11.16) 

where 

a^  +4a_-4a 

4'  *  — - - — 5-  (11.17) 

4(a^a2-a^) 

and  a^,  a^  and  a  are  defined  in  (6.16).  We  note  that  they 
depend  explicitly  on  X^  . 

In  order  to  simplify  our  calculations  of  the  asymptotic 
expressions  for  a^^,  ag  and  a  ,  we  will  restrict  them  to  the 
case  when  X^  *  1  •  Then  introducing  (8.2),  (8.5)j^  2»  and 

(8.8)  into  (6.16),  and  taking  X^  “  1  ,  we  obtain 


47. 


=  4{l  +  +  ^q)  0(n^)}  > 

ag  =  4{1  +  n^(^+  j?) 

a  =  2{l  +  ^cn^  +  n\||c  + JAq)  +  0(n^)}  , 


where 


A.  =  A 


(0) 


and 


^3=1 


W 


(0) 


C  = 


w(o)+wCo!) 
1  2 


With  (11.18)  we  obtain  from  (11.17) 

4-  =  +  *  0(n^)}. 


Then,  introducing  (11.20)  into  (11.16),  we  have 


"IT^o  *  0(n^)} 


(11.18) 


(11.19) 


(11.20) 


(11.21) 


Calculation  of  G[u]  when 


We  now  introduce  into  the  expression  (7.26)  for  G[u]  , 
the  expressions  (11.9),  (11.11),  (11.13)  and  (11.21)  for 
G^,  Gg,  and  G^  and  take  A^  *  1.  With  (8.7)j^  and  (8.2) 
we  thus  obtain 


G[u]  *  +  h^("Tp+  5  ^0  *  ”3~^  “  "9~^  ■'■OCn  )  > 


(11.22) 


where 


48. 


K  = 


=  iwf°)  ♦  =1  • 


(11.23) 


It  is  instructive  to  make  the  substitutions 


£2  =  ,  6 


v,(o) 

16  ”2 


2W 


Co) 


(11.24) 


in  (11.22),  which  may  then  be  rewritten  as 

—  —  2  08  4 

G[u]  -  -  [i-n^(^-6)  +  ocn^l  . 

zil  ^ 


(11.25) 


If  the  material  is  neo-Hookean,  6*0  and  equation  (11.25) 
becomes 


G[ir] 


ill 


(11.26) 


We  may  compare  this  result  with  that  obtained  in  a  previous 
paper  [  2 ]  in  which  an  analysis  similar  to  that  in  the  present 
paper  was  carried  out  for  a  neo-Hookean  material  with  =  j 
and  arbitrary  .  It  was  found  (see  (7.19)  in  [2])  that 

G[ul  -  -773  [1  -4(16 +-^)Ti^  ♦0(Ti'*)I  .  (11.27) 

^^2  ^  2+x3 

By  taking  *  1  in  (11.27)  and  ^  in  (11.26)  we  see 

that  agreement  is  obtained. 


49. 


Appendix  A 


In  this  section  we  shall  prove  the  result  expressed  by 
equation  (4.7). 

We  introduce  the  notation  (c£.(3.3)j^) 

^  (X^-Xg)!  (Al) 

and  the  operator  k[u,u]  defined  by 

<[eu+e^u]  =  e^ic[u]  +  2e^ic[u,u]  +  e^ic[u]  .  (A2) 

With  (2.10)2»  we  obtain 

-N  _  ^  ^  A  _  _ 

•  “a.eVe 

With  this  notation  and  equations  (4.3)  and  (3.3)2  3' 
obtain  from  (2.12),  by  neglecting  terras  of  higher  degree  than 
the  fourth  in  e  , 


2eX  (l-X^)u  +e^{tc[u]  +  2k}  +  2e^ic[u,u] 

1  1,1  ~  -  ~ 

♦  e‘*{ic[u]  +  (x2  +  2x2)e2}  , 


2e2EA  +  4e3EX  (l-x2)G 

c,  ^  X  f  X 

♦  e^E{2k[u]  +  4A, (1-X2)^  +  aI)  . 

X  1^11 


Again  to  order  ,  we  obtain  from  (A4) 


so. 


4E^X^Cl-X^)G^^^{ic[u]+2k} 

+  e^{  (k[u]  +8X^(1-X^)u^  ^k[u,u]  +4K[u]k+4^}  , 

8ehl(l-X^)^(u^^^)^ 

+  12e‘*x2(l-X^)2(G3^^^)^  {ic[u]  +  2k}  , 

16eSJ(l-x2)‘'(G^^^)‘*  ,  (AS) 

4e3A2X^(l-X^)u3^  +  2e‘*E{A2ic [u] +  2A2k 

+  4X2(1-X2)2(a3^^^)2j  ^ 

4£^A^E^  . 

Also,  from  (2.11)  and  (2.16)  we  obtain  to  order  , 

wC3)  =  k^i3  ^  =  kj^i**  ,  (A6) 

where 

>=31  ■  T*?f*112*2X3'»i22*X>222)  .  CA7) 

’‘k  "  H^’*llll**'3’*1112*®^3’*1122**^3*1222*’^3’*2222^  ' 

We  now  substitute  from  (4.3)  in  (2.22),  and  use  the  relations 
(3.1),  (3.2),  (3.3),  (3 .21)  j^,  (A2) ,  (A5) ,  to  obtain,  to  order  e*^, 

G[u]  *  G[eu+e^u] 

*  e^GgLu]  +  2t^t^  II  (gp^+gp^  +  eg^^^)d53^dC2  .  (A8) 


51. 


where 

=  2{k^K[u,u]  +  , 

*  X^(l-x2)[2{2x2w2+(x2-x2)k2+A2k2i}Ea^^^ 

*  AX^Cl-X^)^k2(u^^3^)^]  , 

gf**^  =  k^{<[u]  +  (x2+2x2)e2} 

+  X2W2{2k[u]E  ♦  4X^(1-X^)Eu3^^3^+A^^} 

+  k2{JCic[u])2  +  2X3^(1-X^)G^  .Lic[u.u]+kic[u]+i^} 

+  k2i{4X^(l-x2)2(Ui^l)^E  +  A2ic[u]E  +  2A2ICE} 

+  2X^W22A2E^  +  12x2(l-x2)2k^(G^^^)2{K[u]  +  2k} 

.  8x2(l.x2)2A^k3^(a^^,)2E 

+  16xJ(l-x2)‘*ki^(G3^^^)‘*  . 


CA9) 


52. 


Appendix  B 

We  substitute  in  (5.4)  the  expressions  for  u^,  u^  and  p 
given  in  (3.14)  and  obtain,  with  (3.16)  and  (3.24), 

-Xj^(X^-l){2X^W2-(x|-X^)k2.A2kjj,)E 
Fgi  *  '*'  ^^2,1  • 

fl2  ■  k^C5i,2*kil2_i)-(-l)"«i2“”“«l  *  X“l.2  • 

^22  •  ^l(^2,2-^^l.l^-^2^P  ■  » 

-X^(x2.1){2x2w2-(x2.x2)k^.A2k2^}u^^^ 

-f33-(-l)%33COs2nC3^  , 

where  the  f's  and  cji's  are  functions  of  only  defined  by 

fll  *  (4XX2n2)-X{2k^3’u’ 

+x^(x^-i)k2[3Cx^+i)n%’^+u"^+x^n'*u^+2x^n^(uu’) '] 
+24xJ(x2-l)3k3n2u'2}  , 

fgg  “  -(2x2x2n2)-X{k^8'u'+2x2(x2-l)k2x2n2u’2}  ^  ^32) 

f33  =  -C2x2£22)-l([x2w2-^(x2-x2)k24A2k2i] 

X[(X^+1)  n2u'2+u'»2^^2fjl*u2+2x2n2(;uu’)  ’] 

+  2X^CX^-l)^Ik:23^-6CX2-X3)k3+2A2k3i]n^u'^}  , 


53. 


and 

^11  "  (4AX2“^)'^^2k^6'u'+x2(x2.i)]c^[3(x2+l)£j2u.2 

-U"^-X^a^V^-2X^Q^a]  *  24A!j(A^-l)^k^J22u'^}  , 

*^22  “  ^22  ’ 

4)33  =  -(i2X^n^r^{[X^^V^-\ixl-X^^n^*^h^k^^]  (B3) 

x[  (A^+l)n^u'  2-u"^-A^nS^-2A^n^a] 

+  2A^(A2-l)2[k2^-6(A2-A2)k3+2A2k3^]n2u’^}  , 

<^21  =  (2A^A2n^)'^{k^e’u"+A^(A2-l)k2A^n^U' (U’Vfi^U)}  , 

^12  “  C2AA2n)"^{k^e’u+A^(A2-l)k2U’ (U"+A^£2^U)}  . 


ife  introduce  the  notation 


4a^a 

<l>’  =  -  V  ■—  t 


6aJ(A+1)\3 
Yl -  ^21  ’  ®  “  - iq - 


and  (c£. (3.18) ) 


A  =  A 
A^ 


(A'^-Dk, 

i 

XX^k^ 


Then,  we  obtain  from  (B3)j^  2  4  5  ’  ^^^^2  C3.17) 


11 


22 


■  ,  {28’U'+A[3CA2+I)n2u’2.u"2 

41X2^ 

-A2n‘*u2-2A2jj2a]  +  4  (A-l)  3Bn2u '  2}  , 

- C6'u’  +  2Aa2q2u’2) 


(B4) 


(B5) 


(B6) 


> 


54. 


4>i2  =  {3'U+au' (u"+x2n^u)}  , 

(|)'  =  [u"^-{i+(A^-i)A}n^u’^+ AA^o^Cu'^+n^u^)]  ’  . 

We  note  from  (3.17),  (3.19)  and  (3.25),  with  the  notation 
(B5) ,  the  relations 

<x  =  UU"-U'^  ,  a'  =  UU'" -U'U"  , 
a"  =  n^{A2+l+(A^-l)A}(a+U'2)-(U"^+X^n**U^)  , 

(B7) 

6  *  U"-{l+(X^-l)A}ft^U  ,  e'=U"’ -{l+(x2-l)A}fl^u'  , 

B"  *  X^J2^(U"-n^U)  . 

With  these  relations,  equations  (B6)2  3  yield 

♦22  '  2  [U"^+{x2-l+(3X^+l)A}n2u’2.x2nV]  '  , 

4x^X2^'^ 

*^12  "  ITTl?  Ul+A)(U’u"'+x2RW')-X^n‘*u2  (B8) 

+  AU"2-(l-A)n^U'2}  . 

We  now  substitute  from  (B6)j^  and  (88)2  in  the  expression 
(6.9)  for  q  and  use  the  relations  (B7)  to  obtain 

q  =  [{x2.3+(X^+3)A}a2a-(l-4A)u"2 

+  {x2+3-(x2+9)A-8(X-l)3B}n2u’2 

-(3-2A)x2n‘*u2-2(l-A)U’u"’ +2(1+A)x2n2uu"]  .  (B9) 


From  (3.19)  and  (B5)  it  follows  that 


55. 


2(U'U'")’  =  [(U")^+{X^+l+(x2-l)A}n2u'2-A2j2'*u2] (BIO) 


Then,  from  (B9) ,  (BIO),  (B8)j^,  (B6)^  and  (B4)j^  we  obtain  the 
following  expression  for  4  defined  by  (6,9): 


i  =  ^X"^[3{(X^-1)  +  (X^+l)A}Q^a 

-{5(X^+1)A-(X^-1)A^+8(X-1)^B}S^U'^ 

+  5A(U"^+X^n**U^]  .  (Bll) 

From  (3.20),  (B5)  and  (B7)  we  obtain 
a  =  -(U'2+X^fl%^)  , 

a'  =  {3A2+l+(x2-l)A}n2uu’  ,  on  Cg »  .  (B12) 

B’  =  . 


With  these  relations  and  (3.20),  we  obtain  from  (B6)  and  (B9) 


q 


■j— i  UU  ,  *^22 (1+A)U  ^  , 

2 

^  [2{2(X^-1)+4A-(X^-1)A^+4(X-1)^B}u'2 
+  X^{4x2-(X^-l)A}n^U^]  . 


(B13) 


Equations  (B12)  and  (B13)  yield  the  following  expressions  for 
4^  and  <I>2  defined  in  (6.13): 


$1  =  ^^X“'^{7X^  +  1+(X^-1)A}UU'  , 

$2  =  2fi^X‘^[{2X^-l+2CX^+l)A-j(X^-l)A^ 

+  2(X-l)^B}U’^-^X^(X^-l)An^U^]  . 

From  (B7)g  and  (3.19)  we  readily  obtain,  with  (B5) ,  the 
identities 

X^n^^U^  =  [U' (u"  +  X^n%)-U{U’"  -[2X^  +  1+(X^-1)A]Q^U'}]  ' 
-2X^fl^(UU')  '-U''^-{X^+1+(X^-1)A}«%'^  , 

X^fi^(UU')  '  +  3'U'  =  {U'  (u"+X^fl^U)}'-U'’^-{l+CX^-l)A}fi^u'^  . 


We  now  substitute  from  (B2)j^  2  3  (6.18)  and  (6.17)2 


and  use  (3 . 20) , 

(B4) 

2,  (B5)  and  (BIS)  to  obtain 

b  =  .-1 

f^2 

[u'’^+fi^U'^{l-^(X^+3)A 

-I2 

+  i(x2-i)^-(X-l)3B}]dC2  , 

'>2--  S 

1^2 

-^2 

-[(X^+l)+(X^-l)c]A+(X^-l)^d}]dC2  . 


where  c  and  d  are  defined  by 


-12  1  1  - 
c  =  ^  (XsWg+yAgkgi)  ■  2^  x^~l  ^  ’ 

3 


(k^  +2A^k,J-2X 


X-X; 


TTX:  ""  7^7773 


3  1 


(X  +  1) 


B  , 


defined  in  (3.2),  k^  in  (A7) ,  and  A^  in 


(B14) 

(BIS) 

(B16) 

(B17) 

(3.3) . 


with  k^  and  k^^^ 


57. 
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